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Abstract 



After making some critical comments about the traditional method of ex- 
tracting the gluon condensate from lattice QCD data, I present an alternative 
analysis. The result is more than a factor of five larger than the phenomeno- 
logical value. Two closely related subjects, the effects of the infrared renor- 
malons on the extraction and the Lepage and Mackenzie improvement on a 
lattice perturbation series, are also discussed. 
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One of the popular approaches to study Quantum Chromo dynamics (QCD) in the strong 
couphng region is the QCD sum rule method initiated by Shifman, Vainshtein, and Zakharov 
(SVZ) more than fifteen years ago The success of the approach has been demonstrated 
by many examples ranging from hadron masses, light-front wavefunctions to decay widths, 
form factors, etc. Central to the QCD sum rules is the concept of vacuum condensates which 
are expectation values of composite operators in the QCD vacuum. Among those, the most 
legend is perhaps the gluon condensate defined through the operator = F°'^Faf3, where 
F"^ is the strength of color gauge fields. 

To understand the sum rule phenomenology at a more fundamental level, one has to cal- 
culate the vacuum condensates directly from the QCD lagrangian. Calculations using lattice 
Monte Carlo started in several groups shortly after the publication of SVZ's paper |^. Before 
I comment on these calculations, it is important to point out that phenomenological conden- 
sates from fitting experimental data are in principle different from theoretical condensates 
that are calculated as matrix elements in the QCD vacuum. The former are extracted with 
the Wilson coefficients computed to a few loops, and thus may contain large uncalculated 
multi-loop contributions from the coefficient functions and higher-power corrections. 

The traditional approach of calculating the gluon condensate (in the quenched approxi- 
mation) goes like this Consider an elementary plaquette on the lattice. Calculate using 
Monte Carlo the trace of the plaquette as a function of the lattice spacing a, or the lattice 
coupling constant /3 = 2Nc/gQ. For small enough a, one has. 

The various terms on the right-hand side have different characteristic f3 (or a) dependences. 
The leading term is a series logarithmic in a (or power-like in The condensate term is 

quartic in a (or exponential in By calculating both the left-hand side and the perturba- 
tion series for a wide range of j3, the condensate can be extracted by fitting the expected f3 
dependence. 

The approach shall work in principle. But it is difficult to implement in practice. Let 
me list a few of the practical problems with the approach: 

• First, a big lattice is required to compute the plaquette average for small a. Physically 
the gluon condensate comes from the effects of long wavelength gluons (somewhere 
from 0.5 to 1 fm). The combination of a small lattice spacing and a reasonable physical 
size requires a big lattice. In ref. it was determined that the asymptotic region, 
where the data behave like a sum of powers plus an exponential in jS, starts from f3 = 
6.58. Using a two-loop relation between (3 and a, a = (/5/66o)^^'^^''°/^i6xp(— /5/126o) 
with Al = 4:A MeV, bo = ll/ievr^ and h = 102/25671^, one finds that the coupling 
corresponds to a lattice spacing 0.055 fm. Thus a lattice with 8 points in spatial 
directions spans a physical dimension of 0.44 fm. For f3 = 7, the lattice size further 
reduces to 0.27 fm. Both lattices seem to be too small to measure non-perturbative 
physics. 

• Second, relative to the leading term, the condensate term contributes little to Eq. (1) 
as a decreases. At a = 0.055 fm, the ratio between the leading and the condensate 
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terms is more than 10^. Thus to determine the condensate with a reasonable error at 
small a, one has to compute the perturbation series to large orders. 

• Finally, the perturbation series in Eq. (1) actually divergences due to the infrared (IR) 
renormalons, that is, the coefficients of the series increase like n\ and with a fixed sign 
[^]. The presence of the renormalons complicates the extraction of the condensate in 
two ways. 1) The perturbation series cannot yield better accuracy beyond a certain 
order. 2) The difference between power and exponential behaviors in (3, which has 
been the basis for fitting the condensate, disappears to a certain degree. 

The physical significance of the gluon condensate in light of the IR renormalons has 
been discussed in the literature for a long time 0-0]. The upshot of these discussions is 
that the gluon condensate is a procedure-dependent concept. [In this sense, the status of 
a theoretical condensate is not much different from a phenomenological one.] When a con- 
densate is calculated in a particular scheme designed to regularize the infrared renormalons, 
the coefficient functions in an operator product expansion must be calculated accordingly. 
A consistent method to regularize the renormalons both in the coefficient functions and 
condensates, generalizing the one-loop discussion in ref. 0, has recently been studied by 
this author p. 

In light of the above observations, I consider an alternative strategy to extract the gluon 
condensate from Eq. (1). Assuming the plaquette perturbation series is an asymptotic 
one, one would expect that the magnitude of its terms initially decreases, then reaches a 
minimum, and finally increases without bound. It is generally believed that the minimal 
term occurs at order n ~ /5, with a magnitude ~ a^. Thus the minimum uncertainty in 
the perturbation series induces an uncertainty in the gluon condensate, which can only 
be eliminated through a regularization of the series. However, if the uncertainty is small 
compared with the condensate itself, as Novikov et al. argued 0, it still makes a good deal 
of physical sense to extract a gluon condensate independent of the regularization scheme. To 
illustrate this more clearly, I schematically show in Fig. 1 the order at which the minimum 
term occurs as a function of f3 (the solid line labelled by M). I also show the order at which 
the size of the perturbative term is roughly that of the gluon condensate (with the solid line 
labelled by C). If the uncertainty on the condensate induced by IR renormalons is small, 
the C line would be significantly below the M line. 

Now suppose the perturbation series has been calculated to some fixed order uq, shown 
by the dashed line in Fig. 1, and call its intersection with curve C, /3max- Clearly, one 
cannot extract the gluon condensate accurately ii (3 > On the other hand, for a small 

/3, the corrections from the lattice discretization and the higher-dimensional condensates in 
the expansion become important. Thus there exists a lower limit on /5, /3min; below which 
the extraction becomes unreliable. The existence of a window in [3 (/?min < /?max) depends 
on hq. So long as a window exists, one shall be able to make an approximate extraction. 
The error on the condensate depends on hq. When no is greater than n(/3inin) on the M 
curve, the error is limited by the renormalon singularity. 
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Pmin Pmax p 

Fig. 1. See the text for the explanation of various curves. 

The perturbation series for the elementary plaquette has initially been calculated to 
three loops by Alles, Campostrini, Feo, and Panagopoulos and recently to eight loops 
by Di Reno, Onofri, and Marchesini using a numerical method With uq = 8, /3max is 
somewhere around 6.5 to 7. Experiences with lattice calculations indicate a /^min around 
5.7. In the following analysis, I choose three lattice couplings: /3 = 5.7, 6.0, and 6.2, 
corresponding to the two-loop lattice spacings a = 0.15, 0.105, 0.084 fm, respectively. 

For the average of the elementary plaquette, I use the following result from the lattice 



Monte Carlo [|,|l|,|2 



— TrP= 

Nr. 



0.4509 (/3 = 5.7) 

0.4058 (/3 = 6.0) 
0.3861 (/3 = 6.2) 



(2) 



Since the numerical errors are small compared with the uncertainty in perturbation series, 
I have ignored them. The perturbation series from Ref. M] is. 



^ Cn 2 1.218 2.960 9.28 34 135 563 2488 
> — = — \ \ \ \ \ \ \ h 



(3) 



where I have also ignored the small Monte Carlo errors. For /?'s under consideration, the 
last term calculated is not yet the minimal term in the series. I assume for the moment the 
error in truncating the series is determined by the last term. 



5^-^= 0.4278 ±0.0022 (/? = 5.7) 

n P 

= 0.3988 ±0.0015 (/5 = 6.0) 

= 0.3818 ±0.0011 (/5 = 6.2) 



(4) 



Then the differences between the Monte Carlo data and the perturbation series are, 0.0222 ± 
0.0022, 0.0075 ± 0.0015, 0.0045 ± 0.0011 for p = 5.7, 6.0, 6.2, respectively. In producing 
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the small differences, the knowledge of four, five, and six-loop terms in the series has been 
essential. Notice that the relative error is larger for larger (5. This reflects the observation 
made early — to extract the condensate at larger (3, one needs higher precision for the series. 

The error estimate in Eq. (4) is probably optimistic. One indication is that terms 
in the perturbation series have the same sign and decrease slowly with increasing n. To 
improve the extraction, one may estimate higher-order terms in the series. In the following, 
I shall consider two such estimates: the Pade approximation and the leading IR renormalon 
approximation. 

In the standard Pade approach, both [3,4] and [4,3] approximations produce ([m,n] is 
the standard notation for the Pade approximation with a m-th order polynomial in the 
numerator and a n-th order polynomial in the denominator). 



oo 



E 0.4322 (/? = 5.7) , 
= 0.4011 (/? = 6.0) , 

= 0.3833 ip = 6.2) . (5) 

The differences between the Monte Carlo data and the Pade approximation are, 0.0187, 
0.0047, 0.0028 for P = 5.7, 6.0, 6.2, respectively. The numbers are about two standard 
deviations away from the eight-loop result, indicating that that the eighth-order term may 
not be in the asymptotic region. 



Using the leading IR renormalon series constructed in Ref. [|10|, I find the minimal term 
for the plaquette series occurs somewhere around the 25th to 30th order for the present /?'s. 
The first few terms beyond the eighth order go like this. 



oo 



^ Cn 1.1 X 10^ 4.8 X 10^ 2.1 X 10^ 9.5 x 10^ 4.2 x 10^ 
>— = ... H \ \ \ h 



„ "■■ ' (3^ ' ^11 ^12 ^13 

1.9 X 10^ 8.9 X 10^ 4.1 X 10^ 1.9 x 10^ 9.4 x 10^ 

+ ^li + ^TB + pTe + pl7 + pTs + - U 

Summing the series in the renormalon approximation to the minimal term yields, 

Y: %= 0.4362 (/? = 5.7) , 

n P 

= 0.4029 (/? = 6.0) , 

= 0.3848 iP = 6.2) . (7) 

The differences between the Monte Carlo data and the perturbation series in the leading- 
renormalon approximation are, 0.0147, 0.0029, 0.0013 for P = 5.7, 6.0, 6.2, respectively. The 
numbers are about three standard deviation away from the sum to eight loops. 

As a final result, I take the average of the two large-order estimates as the central 
value and take their difference as the error estimate. Thus the condensation contributions 
to the plaquette average are 0.0167 ± 0.0040, 0.0038 ± 0.0018, and 0.0021 ± 0.0015 for 
P = 5.7, 6.0, 6.2, respectively. These numbers can be converted to the standard form of the 
gluon condensate in unit of GeV^, 
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(^F2)= 0.18 ±0.04 (/3 = 5.7) , 

TT 

= 0.17 ±0.08 {(3 = 6.0) , 

= 0.23 ±0.17 iP = 6.2) . (8) 

Again, the large error for j3 = 6.2 indicates that the perturbation series must be calculated 
to higher orders in order to extract the condensate at large (3. 

The numbers above are consistent with that extracted in Ref. [Q. They are at least 



a factor of five larger than the phenomenological determination (See Ref. |13] for a dis- 
cussion about the phenomenological condensate). Several factors may contribute to this 
discrepancy: 1) higher-orders in the plaquette series are different from what we expected, 2) 
quenched approximation, 3) underestimation of the lattice spacing, 4) the phenomenological 
condensate is contaminated by higher-order Wilson's coefficients. Clearly, it is interesting 
and also important to resolve the discrepancy in this high-precision confrontation between 
lattice Monte Carlo, QCD perturbation theory, and hadron phenomenology. 

The perturbation series for the elementary plaquette (Eq. @) is the first that has been 
calculated up to more than two-loops in lattice QCD. Consequently, it is interesting to 
explore its convergence properties. In general, because of the large tadpole contributions, 
lattice perturbation series converges very slowly. The plaquette series does show such slow 
convergence due to an abnormal fast increase in the coefficients, as observed by Di Renzo 
et al. If a perturbation series has been calculated to sufficient large orders, the slow 

convergence is not a problem. In fact, as I will argue late that the lattice perturbation series 
has the virtue of producing small renormalon uncertainties. 



Lepage and Mackenzie have recently pointed out that the perturbative expansion 
with a more "physical" coupling accelerates the convergence of a lattice series. Clearly, 
if a series is limited to the first few terms, the acceleration improves its predictive power 
considerably. The plaquette series provides an excellent example for studying the effects of 
acceleration. To show this, let me introduce 

F = ^(l-— TrP)-l . (9) 



The series for F is 



F = 1.2755 a^o ± 6.4920 a% + 42.6276 a^Q + 327.1 ajo 

±2720.1 a% ± 23758.8 a% ± 219899 aJo ± ... , (10) 

where aso is the bare lattice coupling. If truncating to the leading term, the result for F, 
0.1015, is only half of the full series, 0.2087, at /5 = 6.0. The bad approximation by the 
leading term is related to the large coefficient 6.4920 in the second term and, in general, 
the fast increase of the perturbative coefficients. However, If changing to a renormalized 
coupling which coincides with the one-loop MS coupling at scale fx = vr/a, 

a.o = as/(l±3.88 a,) , (11) 

F becomes, 

F = 1.2755 a, ± 1.543 ± 11.44 a'l ± 49.7 + 265.8 ± 1566 ± 9616 aj ± ... . (12) 
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Clearly, the convergence has been improved considerably due to the small expansion coef- 
ficients. The leading term is now 0.1584, a number much closer to the full result. This 
phenomenon has already been studied in Ref. . 



Two comments can be made about the acceleration of the lattice perturbation series, both 
of which are related to the question of precision, important, as we have seen, in extracting the 
gluon condensate. First, in the expansion with a new coupling, the coefficients of the lower- 
order terms in general do not increase at a monotonic rate, i.e., they become less regular. 
Regularity is a necessary (not sufficient, of course) diagnosis for the series in the asymptotic 
region, where one can estimate the error of truncation using the last term included. To 
produce a series with regular lower-order terms, one must use some "natural" couplings 
with well-defined relation to the lattice coupling up to high orders. One possibility is to use 
couplings associated with some "physical" schemes like MS. However, the relation between 
the lattice and MS couplings is difficult to calculate. Only recently, the relation has been 
computed to two-loops by Luscher and Weisz [T^], 



a,o = as-3.88 al + 7.7Ual + .... , (13) 

where is again defined at scale fi = n/a. Using a rationalized version of the above relation, 
aso = as/(l + 3.88as + 7.310a^), I have, 

F = 1.2755 as + 1.543 + 2.128 + 29.20 + 169.0 + 880.5 + 9835 aj + ... . 

(14) 

The coefficient of term is now reduced significantly, compared with that in Eq. (p2|). 



However, the irregularity has been shifted to higher-order terms. This is refiected by the 
fluctuation of the successive terms, which makes error-estimate difficult. 

The second comment is about the intrinsic error of the series due to IR renormalons. 
In new expansions, the series converges faster, however, the minimal term in the series in 
general becomes larger. Take the above series as an example, the smallest term in the series 
is approximately 0.0030. Before the improvement, however, the smallest term in the series 
is less than 10~^. This dependence of the renormalon uncertainty on expansion schemes is 
not well studied and is quite interesting. It may imply the existence of a scheme in which 
the coupling constant approaches to zero and the power series converges inflnitely slowly, 
but, the renormalon ambiguity may be vanishingly small. 

To summarize, I have critically examined the traditional method of extracting the gluon 
condensate from the lattice QCD data. Using the plaquette series calculated recently to 
eight loops, I have made a direct determination of the condensate at low (3 and large lattice 
spacings. The result is at least a factor of flve larger than the phenomenological condensate. 
Although still higher-order terms in the perturbation series must be computed to reduce 
the uncertainty, the intrinsic error due to IR renormalons seems to be small. The plaquette 
series converges faster when expanded with more physical couplings, however, the estimation 
of errors becomes more difficult. 
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